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ABSTRACT 

A method for inference of the primordial power spectrum from the COBE^ 
DMR sky maps is discussed. This approach involves a Fourier decomposition of 
the sky maps in a basis of orthonormal functions on the incompletely sampled sky, 
a detailed mathematical model of all noise properties, and a likelihood analysis 
of the data based on an exact probabilistic model of the measurable quantities, 
valid in the context of Gaussian theories of structure formation. Results from the 
two year GORU-DMR data analysis and the implications for inflationary (flat 
geometry) models dominated by either dark matter or a cosmological constant 
are presented. The cross power spectrum of the 53 and 90 GHz DMR sky maps 
is also discussed. 


1. INTRODUCTION 

The development of inflationary ideas during the 1980s (Blau & Guth 1987) 
induced a decade-long adherence to the cosmological paradigm which posits that 
the universe is spatially flat. Such a picture requires that the present energy den¬ 
sity of the universe is dominated by non-baryonic dark matter or alternatively 
by a non-zero vacuum energy contribution (a cosmological constant term. A). 
The minimal version of the model, which invokes cold dark matter (CDM) as 
the major constituent of the universe, is presently in direct confrontation with 
astronomical observations. An extension of the model, which in addition to CDM 
postulates an admixture of hot dark matter (HDM), enjoys considerable popu¬ 
larity in contemporary cosmological research. Vigorous discussion ensues in the 
literature as to the plausibility of the mixed dark matter (MDM) model as a viable 
cosmology, (see e.g. Schaefer, Schah & Stecker 1989; Davis, Summers, & Schlegel, 
1992; Klypin et ah 1993; Pogosyan & Starobinskii 1993; Ma & Bertschinger 
1994; Primack et ah 1994; Ma, and Klypin et ah in this volume). Cosmological 
constant dominated, spatially flat, cold dark matter models (CDM-A ) can be con¬ 
sidered a natural extension of the CDM scenario. By adjusting the combination 
of a non-zero cosmological constant and the total matter density parameter, Oq 
( which remains poorly constrained by observations), we can maintain the spatial 
flatness of the universe required by the inflationary paradigm, whilst striving to 
improve the agreement of the theoretical model predictions to the observational 
data (Peebles 1984; Ffstathiou, Maddox, & Sutherland 1991; Kofman, Gnedin, 
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Bahcall 1993; Carrol, Press, Turner 1992). The empirical determination of the 
power spectrum of primordial inhomogeneities and assessment of their consistency 
with the predictions of these inflation-related models of the universe are critical 
issues in contemporary cosmology. 

The implications of recent observational advancements in extragalactic as¬ 
tronomy, including measures of the galaxy distribution and bulk flow motions for 
the above models of the universe are discussed elsewhere in this volume. In this 
contribution I focus on some issues related to the current status of the cosmic 
microwave background (CMB) anisotropies. White, Scott, & Silk (1994) have 
recently reviewed both the theoretical and experimental status of the held — the 
interested reader is invited to refer to this paper. In recent years there has been 
rapid growth in the experimental area of medium and small angular scale CMB 
anisotropy, but the held is presently in a transitory stage. Therefore I shall con¬ 
centrate on the COBE DMR results and their relation to Dark Matter in the 
universe. 

The COBE DMR discovery of CMB anisotropy (Smoot et al. 1992, Bennett 
et al. 1992, Wright et al. 1992) has affected cosmology in both ontological and 
practical ways, but its predominant quantitative inhuence has been to provide the 
means for the accurate normalisation of theories of large scale structure formation. 
In the following I will discuss the issues of power spectrum estimation from the 
COBE DMR sky maps, the related issue of model normalisation and implications 
for the inhationary models of structure formation. 

2. CMB ANISOTROPY POWER SPECTRUM ESTIMATION 
2.1 COBE DMR Data 

All the results discussed were derived by analysis of the COBE DMR two 
year sky maps. A summary of two years of observations has been presented 
in Bennett et al. (1994). The COBE DMR instrument has observed the mi¬ 
crowave sky at three frequencies (31.5, 53, and 90 GHz) with pairs of ~ 7° 
(EWHM) antennas separated by 60° on the sky. The temperature difference 
measurements were used to construct whole sky anisotropy maps (Janssen & 
Gulkis 1992) binned into 6144 equal area (~ 2?6 X 2?6) pixels. The 31 GHz 
maps were not used in the analysis, since they are relatively noisy and contain 
a more signihcant galactic signal contribution. The A and B channels of the 
sky maps are coadded using inverse noise variance weighting. The sky maps 
have been generated in both galactic and ecliptic coordinate frames, involving 
separate binning of the underlying time-ordered series of sky temperature mea¬ 
surements into 6144 sky map pixels (of equal area). In this analysis, both the 
galactic and ecliptic coordinate frame data sets are used as a check on the extent 
to which the coordinate dependent noise binning can affect the inferred normal¬ 
isation. Strong galactic plane emission was excised from the maps by removing 
all pixels within \b\ = 20° of the galactic equator. Eaint, high latitude galac¬ 
tic emission was not removed from the data (since no fully supported models 
or dehnitive measurements thereof exist). As a simple measure of the extent to 
which such foreground emission could affect the inferred normalisation, the power 
spectrum analysis is conducted with the lowest order (7 = 2) anisotropy mode, 
the quadrupole, either included or excluded. No smoothing or hltering of the sky 
maps is applied prior to the implementation of the power spectrum inference. 



2.2 Model Power Spectra of CMB Anisotropy 

The cosmological models are specified to within an arbitrary amplitude of 
the perturbations as follows: 

1) the global geometry is flat with Oq + = 1, with the bulk mass den¬ 

sity provided by either CDM or MDM; for the MDM models (with only A = 0 
considered) the hot dark matter is introduced in the form of either one or two 
(equal mass) families of massive neutrinos, with the contributed fraction of crit¬ 
ical density taken as = 0.15, 0.2, 0.25 and 0.3 for one massive flavour, and 
Oj, = 0.2, 0.3, otherwise; 

2) the value of the Hubble constant, Hq = 100 h km s“^ Mpc“^, is sampled at 
h = 0.5, and 0.8, and, in agreement with the Big-Bang nucleosynthesis arguments, 
the baryon abundance obeys the relation = 0.013(Reeves 1994); 

3) random-phase, Gaussian, scalar primordial curvature perturbations (no grav¬ 
ity waves) are assumed with the inflationary Harrison-Zel’dovich spectrum corre¬ 
sponding to an adiabatic density perturbation distribution, P{k) oc k. 

The angular distribution of CMB temperature anisotropy induced by such 
curvature perturbations (Sachs & Wolfe 1967, SW) is easily described in Fourier 
language, using the spherical harmonic decomposition of the random held of tem¬ 
perature huctuations — 6TlT[6^(f>) = (f)) (e.g. Peebles 1981). In¬ 

dividual spherical harmonic coefficients aim are Gaussian-distributed in the the¬ 
oretical ensemble of initial conditions for structure formation. The variances of 
the probability distributions of individual modes, (|a^mP), are uniquely expressed 
as integrals over the power spectrum, and depend only on £ due to the statis¬ 
tical isotropy of the CMB temperature held. The CMB anisotropy multipole 
coefficients and the matter perturbation transfer functions for all models were 
evaluated using the Boltzman equation code of Stompor (1994) by solving the 
propagation equations up until the redshift z = 0. 

Over the low-f range of CMB multipoles probed by COBE-DMK the theo¬ 
retical spectra are indistinguishable for the CDM and MDM models with equiva¬ 
lent h and O^. Thus, the power spectrum amplitude derived from the data applies 
equally to both CDM and MDM models. The low-f shape of the CDM/MDM 
spectrum is a little steeper than the underlying, inhationary n = 1 spectrum 
(Bond 1993; Bunn, Scott, & White 1994; Gorski, Stompor, & Banday 1995). The 
steepening is caused by three effects: 1) the anisotropy of the photon distribu¬ 
tion function at the moment of decoupling, 2) the high-redshift integrated SW 
effect generated by the change in the growth rate of matter perturbations at the 
post-recombination epoch due to the non-negligible contribution of radiation (i.e. 
photons and massless neutrinos) to the total energy density, 3) the contribution 
of the hrst adiabatic/Doppler peak to the low £ multipoles. These contributions, 
though small compared to the overall amplitude of the usual curvature pertur¬ 
bation driven SW effect, should be accounted for in an accurate power spectrum 
normalisation based on the C'ORif-DMR anisotropy data. The CDM-A radiation 
power spectra at low-f ( A 10), with a characteristic enhancement of the lowest 
order modes, are determined by both the usual SW and the cosmological constant 
induced, integrated SW (ISW) effects (Kofman & Starobinskii 1985, Gorski, Silk, 
& Vittorio 1992), and depend very weakly on the Hubble constant. 

We express the power spectrum amplitude in terms of Q^ms-PS^ the value 
of the exactly computed quadrupole, 02 , multiplied by (5/47r)^/^, which depends 
non-trivially on cosmological parameters. This is a generalisation of the normal¬ 
isation introduced in Smoot et ah (1992) for pure power law model spectra. For 



Oo ^ 0.5, the ISW contributions boost the quadrupole over its pure SW value. 
For Oo ~ 0.5 — models with a negligible cosmological constant induced ISW ef¬ 
fect — the resulting exact quadrupole is lower than its pure SW counterpart. The 
steepest spectra arise in models near critical matter density (Oq — 0.8). 

2.3 Power Spectrum Parameter Inference 

The power spectrum estimation method described in Gorski (1994) was ap¬ 
plied. Orthogonal basis functions are constructed from spherical harmonics with 
the exact inclusion of both pixelisation effects and galactic plane excision, which 
is coordinate system specihc, and leaves 4016 and 4038 pixels in the galactic 
and ecliptic sky maps, respectively. [Note that unlike the Karhunen-Loeve, or 
signal-to-noise eigenmode approach (Bond 1993, Bunn and Sugiyama 1994) this 
approach utilizes mathematical constructs which are entirely independent of the 
properties of the theory which one wishes to test with the data.] The resulting 
hnite Fourier series decomposition of the signal over the spectral range f G [2, 30] 
is then performed in this basis, and the Fourier mode amplitudes of the CMB 
anisotropy are used in the construction of an exact Gaussian likelihood function 
for a given theory. As explained in Gorski (1994), this method is unique in allow¬ 
ing the algebraic elimination of the non-cosmological, f = 0, 1, modes (without 
affecting higher-f modes at all) from the power spectrum inference procedure. 

The correlation matrices for both theoretical signal and map noise are eval¬ 
uated using the orthogonal mode expansion, i.e. explicitly reflecting the mode¬ 
coupling due to the galaxy cut. The noise is considered to be sufficiently described 
by a pixel-to-pixel uncorrelated random process (the 60° correlations have a neg¬ 
ligible effect on the inference from the DMR data, see Lineweaver et ah 1994). 

Monte Carlo testing has demonstrated that this method of power spectrum 
inference is statistically unbiased (Fig. 1). 

2.4 Results of Qrms-PS Fitting and Large Scale Structure Statistics 

The actual two parameter contour plot of the likelihood derived in the pure 
power law power spectrum ht to the two year data is shown in Gorski et ah 
(1994), and in Bennett et ah in this volume. A detailed presentation of the 
results for CDM, MDM, and CDM-A model htting to the COBE DMR data can 
be found in Gorski, Stompor, & Banday (1995), and Stompor, Gorski, & Banday 
(1995). For comparison one may refer to Bunn, Scott, & White(1994), and Bunn 
& Sugiyama (1994). A typical one parameter likelihood ht of a hat, CDM/MDM 
model to the two year GORA-DMR data yields a ~ 13(7 signihcant determination 
of Qrms-PS ~ 20 yK. For the CDM-A models the best ht amplitude grows with 
decreasing Oq up to ~ 26 yK. Systematic shifts in the central value of the ht are 
observed due to: 1) differences in the noise pixelisation in the galactic and ecliptic 
coordinate frames, which result in a ~ 0.8 yK difference between the inferred 
normalisation amplitudes, with higher values obtained from the ecliptic maps; 2) 
exclusion of the quadrupole, which produces a ~ +0.4 yK variation in the htted 
amplitude; 3) the uncertainty in the values of h and 0;, (rehected by differences 
in spectral shape over the Grange accessible to DMR) causes an additional small 
spread, ±0.25 /rK, of the htted amplitudes. 

A convenient summation of the proposed overall normalisation for the CDM 
and MDM models is Qrms-PS = (20 ± 1.52)(±0.4 ± 0.2 ± 0.25)/rK, and for the 
CDM-A models Qrms-PS ^ {[20 + 9.1h°-® exp(-02/i0.3/o.057)] x (1 ± 0.077) ± 0.4 
±0.2 ± 0.25}/rK. The error ranges represent the statistical error and uncertainties 
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Figure 1: Results from Monte Carlo simulations of the power spectrum inference 
procedure. Each of the 2000 points shows the locus of the maximum likelihood 
parameters derived from a simulated CMB sky map for an input model with n = 1, 
and Qrms-PS = 20 fiK and noise appropriate to the DMR 53 GHz receivers. The 
hlled circle represents the average of the output parameters and demonstrates that 
the inference method is statistically unbiased. The solid line is the approximate 
68% probability density contour derived from the simulated distribution. 


associated with effects 1 through 3 above, respectively. It will be noted that the 
statistical error on the inferred normalisation is considerably larger than the other 
uncertainties. 

Table 1 contains a representative selection of models and the rms values 
of several large-scale structure statistics computed from the matter perturbation 
spectra according to the usual prescriptions given in the footnotes. 

The cold dark matter theory with a standard choice of cosmological pa¬ 
rameters requires a very high normalisation in order to ht the CMB anisotropy 
distribution. Analysis of the hrst year of COBE-DMK data had already sug¬ 
gested that (cr8)mass ~ 1 (Wright et ah 1992, Efstathiou, Bond, & White 1992), 
and this value increases to ~ 1.4 with two years of data and an improved analysis 
technique. Although this normalisation allows the theory to predict POTENT- 
scale velocities (Dekel 1994) with a comfortably high amplitude, it is still not 
high enough to explain the Lauer & Postman (1994) result, and it also results in 
a signihcant overproduction of density perturbations on scales of % 20/i“^Mpc. 
CDM has often been criticised for its poor match to both galaxy and cluster dis¬ 
tributions. Mixed dark matter models manage to circumvent, to a certain degree. 





these same problems by construction — massive neutrinos partially damp the 
perturbations at those length-scales where CDM looks problematic. Among the 
MDM models those with two species of massive neutrino seem to meet the obser¬ 
vational constraints more comfortably (see also Primack et al. 1994). The larger 
free-streaming radius allows for the suppression of the perturbation amplitude on 
larger scales than in models with one massive flavour. This is reflected in the 
decrease of the predicted {(Tg,)mass values. Nevertheless, the proponents of Dark 
Matter cosmology will have to address the viability of the models viz. the con¬ 
tradictory requirements that there should be a bias between the galaxy and mass 
distribution, as suggested by the galaxy cluster abundance argument (White, Ef- 
stathiou, & Frenk 1993) and the galaxy pair-wise velocities, and there should not 
be one, as suggested by the COBE DMR normalisation. 

The low-Oo CDM-A models can not be rejected solely on the basis of the 
COBE-DMK data (Oq A 0.15, 95% conhdence), although the most likely Oq value 
is strikingly close to unity. Conversely, observations of the matter distribution do 
require low values of the total density, with floh^ ~ 0.1. being the favoured value. 
On larger scales, these models predict bulk flows in reasonable agreement with 
POTENT, but dramatically smaller than determined by the Lauer & Postman 
analysis. Therefore, if the latter observation is conhrmed all CDM-A models 
(together with critical density CDM or MDM models) will be found wanting. 

At present, although the position of the CDM-A models is more comfort¬ 
able than that of critical density models, it is not free from potentially fatal 
flaws. Unfortunately, CMB anisotropies on a one degree scale do not offer any 
serious prospects for distinguishing between cosmological constant models and 
other viable scenarios (Stompor & Gorski, 1994; Bond et al., 1994). Hence, deh- 
nite conclusions will have to await more reliable observational data, particularly 
large-scale, deep galaxy surveys capable of unraveling the shape of the galaxy 
power spectrum down to A; ~ O.OlhMpc"^, or a comoving length-scales up to 
~ 600 /i“^Mpc. 


3. CROSS POWER SPECTRUM 

Visualising the faint CMB structure present in the rather noisy DMR maps 
is a difficult task. A straightforward and relatively simple procedure is the compu¬ 
tation of a ‘cross’ power spectrum between two sky maps. Such a construction can 
take full advantage of the multiple channel data acquisition by the DMR instru¬ 
ment. Since the noise properties of the different channels are not correlated, there 
is no noise bias in the amplitude of the quadratic power spectrum coefficients. 

The cross power spectrum presented here is derived by the direct product 
of the individual Fourier coefficients from the 53 and 90 GHz sky maps. Monte 
Carlo simulations are performed to obtain comparison statistical distributions for 
the CDM (MDM) model, since analytic solutions are no longer easily computed. 

Figure 2 displays both the sum (A+B) and difference (A-B) map cross power 
spectra, and clearly demonstrates a signihcant detection of CMB anisotropy by 
the COBE-DMK instrument. The signal-to-noise ratio in the cross-combined 53 
and 90 GHz two year data exceeds unity for £ % 15. The cross spectra derived in 
galactic and ecliptic coordinates appear to be in excellent agreement. This is not a 
trivial fact — the basis functions with respect to which the Fourier decomposition 
was performed are coordinate frame dependent due to the different geometry of 
the applied Galactic cut. Thus Fig. 2 is the best demonstration to date of the 
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Figure 2: The cross power spectra for the DMR 53 and 90 GHz two year sum and 
difference sky maps derived in both galactic and ecliptic coordinate frames. The 
results of 40,000 Monte Carlo simulations of CDM skies normalised to Qrms-PS = 
19.5 fiK and with noise appropriate to the two year sky maps are shown. The 
heavy, solid line shows the mode of the power distribution at each £. Note that 
this is not a ht to the plotted cross-power spectrum. The grey bands represent 
the conhdence regions for the distribution of power at each value of £ as follows: 
light grey — 68% , medium — 95%, and dark grey — 99.8% . 


rotational invariance of the angular spectrum of COBE DMR measured CMB 
anisotropy. 

The overplotted CDM power spectrum distribution (identical to MDM over 
this range of £) appears to be a reasonable statistical descriptor for the mea¬ 
sured CMB anisotropy pattern. The low COBE DMR two year quadrupole am¬ 
plitude (Bennett et ah 1994) is clearly evident, despite the fact that we have 
not attempted to model and correct for residual high latitude Galactic emission. 
Although such modelling would result in a lower estimated cosmological CMB 
quadrupole amplitude, it would still fall within the 95% probability zone under 
the distribution derived from the COBE DMR normalised CDM/MDM model. 
An analogous plot for CDM-A models would only show differences at low £ for 
low values of Oq! the average power in this range would tend to be higher than 
in Fig. 2. Coupled with the relatively low measured value of the quadrupole, the 
relative likelihood of high-A models is lower than that for the matter dominated 
models. 






5. SUMMARY 


The improved quality of the two year C'ORU-DMR data combined with 
reliable power spectrum estimation techniques allows the accurate normalisation 
of cosmological theories. Previously, the variance of C'ORU-DMR temperature 
fluctuations on a 10° angular scale was utilised for the normalisation of the power 
spectrum. Subsequent work (Wright et ah 1994, Banday et ah 1994, and these 
proceedings) has demonstrated that this technique can be unreliable without con¬ 
siderable attention. More appropriate methods take full advantage of the mea¬ 
sured CMB anisotropy power distribution on all angular scales accessible to the 
COBE-DyVK instrument, as discussed in this contribution. 

Whilst it would be safer to await the hnal 4-year COBE results before offer¬ 
ing dehnitive statements as to the viability of theoretical models, one should note 
that the CDM normalisation derived from the two year C'ORU-DMR data does 
appear to be irreconcilably high, while the MDM and CDM-A models have little 
room left for adjustments. Interestingly, theoretically favoured globally flat uni¬ 
verse models notwithstanding, it is possible to construct open universe inflationary 
dark matter models in reasonable agreement with the observations (Gorski, Ratra, 
Sugiyama, & Banday 1995, and references therein). 
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Table 1: Inferred cosmological statistics for spatially flat CDM, MDM and CDMA 
models (with Oq + +^5 = 1, where Oq = ^cdm + and = 0.013/i“^), 

with a COBE-DMK normalisation expressed in terms of Q^ms-PS- Ni, is the 
number of massive neutrino species, and rrii, the neutrino mass in eV. The errors, 
including both statistical (fcr) and systematic deviations, are of the order of 11%. 
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ABSTRACT 

We summarize the results of the Cosmic Microwave Background (CMB) 
anisotropy measurements of the COBE-DMR using the hrst two years of data. 
The inference of angular power spectrum parameters from the DMR data is com¬ 
plicated by the need to omit Galaxy contaminated regions of the sky, preventing 
the spherical harmonics from forming an orthogonal basis set over the remain¬ 
ing sky. The COBE-DMR team has implemented three methods for deducing the 
power spectrum parameters. These include: a two-point correlation function anal¬ 
ysis (Bennett et ah 1994); an analysis in terms of orthogonal modes constructed 
on the cut sphere (Gorski et ah 1994); and an analysis of a modihed Hauser- 
Peebles power spectrum (Wright et ah 1994). These methods for deducing the 
power spectrum are discussed and the results are summarized and compared. 

1. INTRODUCTION 

The large angular scale cosmic temperature fluctuations are believed to re¬ 
flect the gravitational potential fluctuations at the epoch of recombination. The 
amplitude, Q^ms-PS^ and power law spectral index, n, of the primordial anisotropy 
are important cosmological parameters, and serve as initial conditions for for N- 
body and hydrodynamic codes that test dark matter models. We have applied 
three methods for determining the amplitude and slope of the primordial power 
spectrum to the two-year COBE-DMR data. The results from each method are 
presented in detail in Table 1. The three methods yield maximum likelihood val¬ 
ues for n between 1.02 and 1.42 with uncertainties between 0.32 and 0.45. In 
no case is a determination inconsistent with the scale-invariant value n = 1, and 
the steeper determinations tend to be driven by the smallness of the quadrupole. 
All methods return a similar ridge of degeneracy between Qrms-PS and n, though 
the slope of the ridge does vary slightly from method to method. The maximum 
likelihood values for the normalization of an n = 1 spectrum range from 18.2 to 
20.4 /iK, with uncertainties between 1.6 and 2.0 fiK. The n values returned by 
each method vary by about 0.5 cr, while the normalizations for a scale-invariant 
spectrum vary by about 1 cr. 

The two-point correlation analysis and the modihed Hauser-Peebles power 
spectrum analysis employ maximum likelihood hts to quadratic statistics con¬ 
structed from products of maps: a combination of the 53 and 90 GHz data. The 

^The National Aeronautics and Space Administration/Goddard Space Flight Center 
(NASA/GSFC) is responsible for the design, development, and operation of the Cosmic Back¬ 
ground Explorer (COBE). Scientific guidance is provided by the COBE Science Working Group. 
GSFC is also responsible for the development of the analysis software and for the production of 
the mission data sets. 



Table 1: Power Spectrum Analysis of Combined 53 and 90 GHz 2-yr DMR Data 


Parameter 

2-pt function 

t/jijn modes 

Gim modes 


Bennett et al. 

Gorski et al. 

Wright et al 


including quadrupole 


n 

l-42iE)f 

1 99+0.43 

— 

n (marginaP) 

1.42 ±0.37 

1.10 ±0.32 

1 oq+0.34 
1.03'_o.39 

Qrms — PS 

14.3t^i 


— 

Qrms — PS\n=l (/^K) 

18.2 ± 1.5 

19.9 ± 1.6 

— 


excluding quadrupole 


n 


1 

l.UZ_Q 59 

— 

n (marginal) 

1.11 ± 0.40 

0.87 ±0.36 

1 9r+0.40 
l.ZO_o,45 

Qrms — PS 

n.Pli 

20.0t(;+ 

— 

Qrms — PS\n=l (/^K) 

18.6 ± 1.6 

20.4 ± 1.7 

19.8 ±2.0 

ai at pivot (/uK)^ 

ar = 9.5 ± 1.0 

og = 8.3 ± 0.7 

— 


“Bennett et al. and Wright et al. define the marginal likelihood on n to be 
L{n) = maxiv((5|n), while Gorski et al. dehne it to be L{n) = L[Q^n)dQ^ 
where Q = Qrms-PS- 

^An alternate power spectrum normalization, expressed in terms of the multipole 
moment that is independent of n. 

techniques make extensive use of Monte Carlo simulations to construct and cali¬ 
brate the derived likelihood functions. The orthogonal mode analysis employs an 
exact maximum likelihood ht to the linear coefficients Ci from the 53 and 90 GHz 
data, which, in effect, places the signal and noise on a similar footing. Monte 
Carlo simulations are only used to verify that the technique is statistically unbi¬ 
ased. Note that the three methods will tend to respond differently to unmodeled 
effects in the data, such as residual Galactic emission, or systematic artifacts, 
though we believe that any such effects are small. While the three methods re¬ 
turn generally consistent results, small discrepancies do exist and the reasons for 
this are not presently understood. Studies of our methodologies are continuing as 
we prepare to analyze the full four-year sky maps. 

2. TWO-POINT CORRELATION FUNCTION 

The two-point cross correlation function at angular separation a, C(a), is 
the average product of all pixel temperatures with a hxed angular separation: 
C(a) = J2i,j WiWj TA^iTs.j / where the sum is restricted to pixel pairs 

(z,j) for which the pixel angular separation resides in the bin dehned by a, TA,i 
and Ts^i are the observed temperature in maps “A” and “B” respectively, after 
monopole and dipole (and quadrupole) subtraction, and Wi is the statistical weight 
of pixel i. The analysis discussed below is based on the 53 X 90 GHz cross 
correlation function. 

We determine the most likely quadrupole normalized amplitude, Q^ms-PS^ 
and spectral index, n, by evaluating the Gaussian approximation to the likelihood 








Figure 1: Likelihood contours as 
a function of n and Qrms-PS for 
the53(A+B)/2 x 90(A+B)/2 cross 
correlation data. The contours cor¬ 
respond to 68%, 95%, and 99.7% 
conhdence regions, corrected for a 
bias in n of +0.11, and in Qrms-PS 
of —1.5 /iK. The solid curves are 
contours with the quadrupole in¬ 
cluded in the analysis; the peak 
of the likelihood is indicated by 
+ . The dashed curves are contours 
with the quadrupole excluded from 
the analysis; the peak of the likeli¬ 
hood is indicated by *. See Table 1 
for a quantitative summary of these 
likelihoods. 


function 

exp[-|(AC)^M-i(AC 

[(27r)™det(M)]i/2 

Here (AC)^ and (AC) are row and column vectors with entries (AC')^ = C(ai) — 
(%(«,)), and M = ((AC)(AC)^) is the covariance matrix of the correlation 
function. The angled brackets denote averages over receiver noise and cosmic 
variance. We estimate the mean correlation function and covariance matrix, as a 
function of Qrms-PS and n, by means of Monte Carlo simulations that account 
for ah important aspects of our data processing. The hnal results are calibrated 
with additional Monte Carlo simulations to account for any biases introduced by 
the Gaussian approximation or other effects. Figure 1 (after Figure 4 of Bennett 
et al.) shows the resulting likelihood contours as a function of Qrms-PS and n for 
the analyses both including and excluding the quadrupole moment. The contours 
correspond to 68%, 95%, and 99.7% conhdence regions. 

3. ORTHOGONAL MODF DFCOMPOSITION 



Gorski et al. explicitly construct linear combinations of spherical harmonics 
to form orthonormal functions on the cut sphere, following a procedure developed 
by Gorski (1994). The 53 and 90 GHz data are decomposed with respect to 
these modes: T[p) = Y,i where T[p) is the temperature in map pixel p, 

and Ci is the coefficient of the ith mode ffi. For Gaussian power-law models of 
anisotropy, the Ci coefficients have a multivariate Gaussian distribution, thus an 
exact likelihood function for the parameters Qrms-PS and n may be constructed 


_ exp[-|c^ • (CcMB + Cn) ^ • c] 
[(27r)“det(CcMB + CN)]l/2 


where c is the vector of Ci coefficients obtained from decomposing the 53 and 90 
GHz data, and (CcMB + Cn) is the covariance matrix of the spectra c, which 







Figure 2: Likelihood contours as a 
function of n and Qrms-PS for the 
53 and 90 GHz orthogonal mode de¬ 
composition. The contours corre¬ 
spond to 68%, 95%, and 99.7% con- 
hdence regions. The solid curves 
are contours with the quadrupole 
included in the analysis; the peak 
of the likelihood is indicated by 
+ . The dashed curves are contours 
with the quadrupole excluded from 
the analysis; the peak of the likeli¬ 
hood is indicated by *. See Table 1 
for a quantitative summary of these 
likelihoods. 


depends on Qrms-PS and n. See Gorski and Gorski et al. for details. Monte Carlo 
simulations with controlled inputs indicate that the parameters inferred with the 
above likelihood function are not statistically biased with respect to the input. 
Figure 2 (after Figure 2 of Gorski et al.) shows the resulting likelihood contours 
as a function of Qrms-PS and n for the analyses both including and excluding the 
quadrupole moment. The contours correspond to 68%, 95%, and 99.7% conhdence 
regions. 

4. MODIFIFD HAUSFR-PFFBLFS POWFR SPFCTRUM 

Wright et al. solve for the angular power spectrum of the DMR two year 
data by modifying and applying the technique described by Peebles (1973) and 
Hauser & Peebles (1973) for data on a cut sphere. The input data for this analysis 
is a “cross” power spectrum which involves the evaluation of terms of the form 
{GimTA) • {GimTB) where Ta and Tg are two independent DMR maps, Gim is a 
modihed spherical harmonic which, by construction, is orthogonal to the monopole 
and dipole (and, optionally, the quadrupole) on the cut sky, and the angle brackets 
denote an inner product on the cut sky. See Wright et al. for more details. They 
have evaluated a Gaussian likelihood, similar to eq. 1, over a variety of spectral 
ranges: 3 < 7 < 19, 3 < 7 < 30, and 2 < 7 < 30, and for a variety of data 
combinations. After calibrating their results with Monte Carlo simulations they 
quote most likely spectral indices in the range 1.02 to 1.41 depending on the choice 
of data, the spectral range included, and the type of modes used. Their results, 
together with the other two methods are summarized in Table 1. 
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ABSTRACT 


We summarize the recent results on the observed COBE-DyiB two year 
sky rms temperature fluctuations. A ‘cross-RMS’ statistic is used to infer the 
Qrms-PS normalization for a scale-invariant (n = 1) spectral model. The method 
is extended to the normalization of other cosmological power spectra. 


1. INTRODUCTION 


In principle, the observed sky rms on a given angular scale is a convenient 
number to use for the normalization of a particular cosmological model. Wright et 
al. (1994a) used the sky rms temperature fluctuations from the hrst year COBE- 
DMR sky maps smoothed to an approximate FWHM of 10° to determine the 
effective normalization Qrms-PS for the scale-invariant Harrison-Zel’dovich power 
spectrum, P{k) oc A;”, where n = 1 and k is the comoving wavenumber. We 
have updated and extended the analysis of Wright et al. (1994a) using the two 
year COBE-DMK data to infer the normalization for an n = 1 power law model, 
together with a number of other cosmological power spectra. 


2. ANALYSIS 


Banday et al. (1994) have dehned a cross-RMS, C)_rm 5 , between two maps 
a and b as 


= E T‘T‘ y / u-*, (1) 

i i 

where the sums are over ah pixels i surviving the galactic cut (|b| > 20°), and Wi 
is the weight assigned to that pixel. Unit weighting is adopted here. The best ht 
monopole, dipole and, where appropriate, quadrupole are removed from the sky 
maps prior to the computation of the (E)rms- 

A Monte Carlo approach was adopted to generate the (E)rms distributions 
for a grid of Qrms-PS values (with 2500 simulations used for each value of 
Qrms-Ps)- Bach simulation combines a realization of sky anisotropy hltered 
through the COBE-DMK beam (Wright et al. 1994a) with noise realizations 
based on the specihc COBE-DMK channels. The Qrms-PS-dependent statistical 



means, variances and covariances of the 7° and 10° (E)rms are derived from these 
Monte Carlo simulations, and used to construct the gaussian approximation to the 
probability distribution of the (E)rms- This, together with the measured (E)rms 
values, dehnes the likelihood function C{Qrms-Ps)- That the (E)rms provides a 
statistically unbiased estimator of the true cosmological Qrms-PS amplitude has 
been verifed with additional simulations used as test input data. 

We refer the reader to Banday et ah (1994) for full details. Summarizing 
the main results therein: 

• The observed (E)rms is frequency independent and thus consistent with a 
cosmic origin for the temperature fluctuations. 

• The (E)rms values derived from different combinations of the 53 and 90 
GHz data are generally in excellent statistical agreement (Table 1), thus 
the inferred Qrms-PS values are reasonably independent of data selection. 

• Qrms-PS derived from the (E)rms at 7° and 10° including the quadrupole is 
biased low by the small observed GOHif-DMR sky quadrupole amplitude. 
This can be corrected for, after which the estimates for Qrms-PS both 
including and excluding the quadrupole are in excellent agreement. 

• the amplitude of primordial inhomogeneity in the context of a Harrison- 
Zel’dovich n = 1 power law spectral model is Qrms-PS ~ 19 /uK (in good 
agreement with values derived from other techniques — see Gorski et ah 
1994, Bennett et ah 1994, Wright et ah 1994b). 


Table 1: Observed (E)rms values derived from possible 53 and 90 GHz 
combinations. 


(d)RMS combination 

including Quadrupole 

excluding Quadrupole 

r (pK) 

10° (gK) 

r (,iK) 

10° (gK) 

53A(g)53 

90A(g)90B 

53A(g)90A 

53B(g)90B 

53A(g)90B 

53B(g)90A 

53(A+B)(g)90A 

53(A+B)(g)90B 

53A(g)90(A+B) 

53B(g)90(A+B) 

53(A+B)(g)90(A+B) 

(53+90)A(g)(53+90)B 

(53A+90B)(g)(53B+90A) 

44.51^7 

28.4t?;^ 

oo 0 + 7 . 0 
OZ/.O Q j 

39.3tt? 

33.5t^;? 
35.7t^;^ 
35.6tt^ 

32.4tJ| 

25.7i3;4 

30.9i2;7 

9Q q+2-3 

32A^f^2 

31.711^0 

3L3lli 

31.2i^;® 

31.7tl:® 

30.8tll 

3i.2il:? 

30.7tl;^ 

29.9il;? 

43.1+5-2 

o.oi-:I 

27.7lli 

30.2tri 

44.0t®:® 

OO 7+10.8 
OO. 1 _10.0 

37.7t^;? 

36. 

31.91L" 

34.411^ 

34.4tt® 

28.8it^ 

30.6tl^ 

24.711:1 

30.3ll:l 

27.7lft 

30.9111 

30.8t;^;E) 

30.611:1 

29.3t^;E( 

30.6t^;E) 

29.3ll:l 

29.9ll-^r 

29.4t+^ 

28.5tl;? 



MODEL NORMALIZATIONS 


The normalization of other models of cosmological anisotropy has proceeded 
by a detailed reworking of the above. We report the inferred Qrms-PS values for 
a number of flat geometry, critical density, scale-invariant CDM models. 

An exact calculation for the power spectrum of CMB anisotropy in 
this context renders an effective spectral slope (in the sense of a power law 
approximation, P{k) oc , used to generate the multipole coefficients solely 
through the Sachs-Wolfe effect) slightly steeper than n = 1 over the angular scales 
probed by the COBE-DMK instrument (f A/ 15). A value in the range 1.1 - 1.15 
would be most appropriate. The steepening is caused by three effects: 

1) the anisotropy of the photon distribution function at the moment of 
decoupling. 

2) the high-redshift integrated Sachs-Wolfe effect generated by the change in 
growth rate of matter perturbations at the post-recombination epoch due to the 
non-negligible contribution of photons to the total energy density. 

3) the contribution of the hrst adiabatic/Doppler peak to the low £ multipoles. 

Table 2 summarizes the normalization amplitudes inferred from the (E)rms 
values derived from the 53 and 90 GHz (A+B) sky maps, excluding the 
quadrupole. The results obtained when the quadrupole is included are in 
excellent agreement with these once the bias introduced by the low observed 
quadrupole is accounted for (see Banday et ah 1994). The steeper power law 
approximation hts render lower amplitudes than the exactly computed CDM 
power spectra, and should not be considered as sufficiently good descriptors of 
this class of cosmological model (see Gorski, Stompor & Banday 1995). The 
Qrms-PS normalization for CDM models is also ~ 19 /rK. 

Table 2: Qrms-PS normalizations for specihc cosmological models. 


Model 

Qrms-PS 

power law, n = 1 
power law, n = 1.1 
power law, n = 1.15 

CDM {h = 0.5, 0;, = 0.01) 
CDM {h = 0.5, Ofc = 0.03) 
CDM \h = 0.5, Ofc = 0.05) 
CDM \h = 0.5, Vtk = 0.10) 
CDM \h = 0.8, Ofc = 0.03) 

IQAp'i 

18.2i?i 

n.pPs 

18.9PP 

18.811-1 
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